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(54) Method and apparatus for matrix reordering 

(57) A matrix reordering method performs reorder- 
ing of elements of a coefficient matrix created based on 
coefficients of linear simultaneous equations whose so- 
lutions are to be produced by parallel processing of 
processors of a computer in accordance with Gaussian 
elimination. Herein, degrees corresponding to numbers 
of non-zero elements are calculated with respect to all 
pivots included in the coefficient matrix. Then, a first piv- 
ot whose degree is under a threshold (mindeg> a) is 
selected from among the pivots of the coefficient matrix, 
while a second pivot whose critical path length is mini- 
mum is also selected from among the pivots of the co- 
efficient matrix. Replacement of elements is performed 
between the first pivot and second pivot to complete re- 
ordering with respect to the first pivot. In addition, non- 
zero elements, which are newly produced by the Gaus- 
sian elimination of the first pivot, are added to the coef- 
ficient matrix. If a degree or a parameter of the first pivot 
is under a threshold (p), reordering is performed on a 
partial matrix whose elements are not eliminated and 
are selected from among the elements of the coefficient 
matrix in accordance with the nested dissection method, 
so that non-zero elements, which are newly produced 
by the Gaussian elimination of the partial matrix, are 
added to the coefficient matrix. Because the critical path 
length can be reduced as compared with conventional 
techniques, it is possible to considerably reduce a total 
processing time of the parallel processing to secure 
high-speed performance of the Gaussian elimination. 
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Description 

BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] This invention relates to methods and appara- 
tuses for matrix reordering that is preferable for use in 
determination of orders of coefficient matrices of linear 
simultaneous equations. This invention also relates to 
methods and apparatuses for electronic circuit simula- 
tion using the matrix reordering. 

Description of the Related Art 

[0002] Conventionally, there are provided system 
simulation methods that use computers to describe el- 
ements constructing the system as linear simultaneous 
equations, which are repeatedly resolved to analyze op- 
erations of the system. Gaussian elimination is frequent- 
ly used in the system simulation methods and is known 
as one of methods for obtaining solutions of the linear 
simultaneous equations. The Gaussian elimination pro- 
duces solutions by an advance elimination process and 
a back-substitution process with respect to a matrix of 
nxn+1 , which is formed by linearly arranging constant 
terms and coefficient matrices representing coefficients 
for n elements of linear simultaneous equations, for ex- 
ample. The advance elimination process is effected to 
transform the matrix of n x n+1 to an upper triangular 
matrix, while the back-substitution process is effected 
to sequentially produce solutions for the transformed 
matrix from its last row. Matrix reordering or matrix or- 
dering is known as the technique for reordering ele- 
ments (particularly, non-zero elements) of the coeffi- 
cient matrices representing coefficients of the simulta- 
neous equations in accordance with prescribed rules. 
Using the matrix reordering, it is possible to reduce a 
number of times for performing calculations such as 
multiplication and division that are needed for producing 
solutions to the linear simultaneous equations. 
[0003] FIG. 1 4 shows a flow of processes of the con- 
ventional matrix reordering for use in the electronic cir- 
cuit simulation system. In general, the conventional ma- 
trix reordering employ the prescribed method called 
"minimum (or least) degree method". This is shown by 
a book entitled "Matrix Calculation Software" that is writ- 
ten and edited by Dr. Tikara Kokuni and is published by 
Maruzen Co., Ltd. in Japan. The flow of FIG. 14 is pro- 
vided for reordering elements of a coefficient matrix of 
n x n . In step S1 01 , a variable is set to '1 '. In step S1 02, 
a pivot having a minimum degree is selected from 
among pivots that are not determined in order within a 
given coefficient matrix, then, a number of the selected 
pivot is set to a variable p. The coefficient matrix con- 
tains a number of pivots, wherein pivot k (where k is an 
integer ranging from t to N) designates a combination 
of elements of row k and elements of column k. In addi- 



tion, the degree is defined as a number of non-zero el- 
ements that are contained in a pivot from which ele- 
ments on a diagonal line are excluded, wherein the non- 
zero elements counted for the degree are arranged in a 
s lower side of the diagonal line in a column direction and 
in a right side of the diagonal line in a row direction. The 
pivot having a minimum degree is selected from among 
pivots whose numbers are counted from i. In step S1 03, 
new pivot i is determined by replacing its original pivot i 
io with the selected pivot p. In step S104, non-zero ele- 
ments that are produced by effecting Gaussian elimina- 
tion on pivot i are added to the coefficient matrix. Re- 
placement between pivots p and i is to interchange all 
elements of row p and all elements of row i and to inter- 
's change all elements of column p and all elements of col- 
umn i. That is, replacement is effected simultaneously 
on rows and columns of pivots. A series of processes 
described above are repeatedly performed until all piv- 
ots are to be completely determined in order, which is 
20 shown by steps S105 and S106 in FIG. 14. According 
to the minimum degree method, pivots each having a 
small number of non-zero elements are reordered to- 
ward an upper side within the coefficient matrix. This al- 
lows reduction of a total number of calculations such as 
25 multiplication and division in the Gaussian elimination. 
[0004] Using the coefficient matrix whose elements 
are reordered in accordance with the conventional ma- 
trix reordering method, parallel processing is executed 
on the Gaussian elimination process by using multiple 
30 processors in parallel. As compared with successive 
processing using a single processor, the parallel 
processing using the multiple processors does not al- 
ways sufficiently increase the processing speed thereof 
Because, in the Gaussian elimination process, an oper- 
as ational process (or calculation) for a certain row is made 
with reference to result of an operational process effect- 
ed on another row. In the parallel processing of the 
Gaussian elimination process, operational processes 
(or calculations) are respectively made by units of rows 
^0 and are performed in parallel by the multiple processors. 
To achieve high-speed performance of the Gaussian 
elimination process, it is necessary to consider several 
points in advance. A first point is to reduce a total 
number of operational processes that are effected in the 
45 parallel processing of the Gaussian elimination process. 
A second point is to reduce a number of operational 
processes that are effected on each of the rows. As for 
a first row, for example, its operational process is made 
with reference to result of an operational process of a 
50 second row, which is also made with reference to result 
of an operational process of a third row. Hence, it is nec- 
essary to count all of the aforementioned operational 
processes that are involved into the operational process 
being effected on each row. Although the aforemen- 
55 tioned points are very important, the conventional matrix 
reordering is made in consideration of only the reduction 
of the total number of operational processes being ef- 
fected in the Gaussian elimination process. Therefore, 
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the conventional matrix reordering do not always con- 
tribute to high-speed performance in the parallel 
processing of the Gaussian elimination process. 

SUMMARY OF TH E INVENTION 5 

[0005] It is an object of the invention to provide a 
method and an apparatus for matrix reordering that is 
improved in high-speed performance, particularly in par- 
allel processing of the Gaussian elimination process us- 10 
ing multiple processors. 

[0006] It is another object of the invention to provide 
a method and an apparatus for electronic circuit simu- 
lation that use the matrix reordering improved in high- 
speed performance. is 
[0007] A matrix reordering method performs reorder- 
ing of elements of a coefficient matrix created based on 
coefficients of linear simultaneous equations whose so- 
lutions are to be produced by parallel processing of 
processors of a computer in accordance with Gaussian 20 
elimination. Herein, degrees corresponding to numbers 
of non-zero elements are calculated with respect to all 
pivots included in the coefficient matrix. Then, a first piv- 
ot whose degree is under a threshold (mindeg+ a) is 
selected from among the pivots of the coefficient matrix, 25 
while a second pivot whose critical path length is mini- 
mum is also selected from among the pivots of the co- 
efficient matrix/Replacement of elements is performed 
between the first pivot and second pivot to complete re- 
ordering with respect to the first pivot. In addition, non- 30 
zero elements, which are newly produced by the Gaus- 
sian elimination of the first pivot, are added to the coef- 
ficient matrix. If a degree or a parameter of the first pivot 
is under a threshold (p), reordering is performed on a 
partial matrix whose elements are not eliminated and 35 
are selected from among the elements of the coefficient 
matrix in accordance with the nested dissection method, 
so that non-zero elements, which are newly produced 
by the Gaussian elimination of the partial matrix, are 
added to the coefficient matrix. 40 
[0008] If the coefficient matrix is given as a non-sym- 
metric matrix, a symmetric matrix is created prior to the 
reordering by transposition. Then, the reordering is per- 
formed on the symmetric matrix. 

[0009] Because the critical path length can be re- 4$ 
duced as compared with the conventional technique, it 
is possible to considerably reduce a total processing 
time of the parallel processing to secure high-speed per- 
formance of the Gaussian elimination. 
[001 0] In the case of the electronic circuit simulation, so 
the coefficient matrix is created using coefficients cor- 
responding to nodes by which electronic components of 
a given electronic circuit are connected. 

BRIEF DESCRIPTION OF THE DRAWINGS 55 

[0011] These and other objects, aspects and embod- 
iments of the present invention will be described in more 
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detail with reference to the following drawing figures, of 
which: 

FIG. 1 is a flowchart showing an electronic circuit 
simulation method using matrix reordering in ac- 
cordance with an embodiment of the invention; 
FIG. 2 is a flowchart showing details of a matrix re- 
ordering process contained in the electronic circuit 
simulation method shown in FIG. 1 ; 
FIG. 3 is a flowchart showing details of selected 
steps in FIG. 2; 

FIG. 4 is a flowchart showing a matrix reordering 
process that is applied to the electronic circuit sim- 
ulation method of FIG. 1 in accordance with another 
embodiment of the invention; 
FIG. 5 shows a part of a program that describes a 
calculation method of a Gaussian elimination proc- 
ess contained in the electronic circuit: simulation 
method of FIG. 1 ; 

FIG. 6 shows an example of a partial matrix whose 
elements are not reordered; 
FIG. 7 shows a task graph for a Gaussian elimina- 
tion process of the matrix of FIG. 6; 
FIG. 8 is a circuit diagram showing an electronic cir- 
cuit in which electronic elements are connected to- 
gether by way of nodes; 

FIG. 9 shows a coefficient matrix whose elements 
are configured and arranged in connection with the 
nodes of the electronic circuit; 
FIG. 1 0 shows a reordered matrix that is created by 
reordering elements of the coefficient matrix of FIG. 
9 in accordance with the matrix reordering process 
of FIG. 2; 

FIG. 11 shows a task graph designating a critical 
path that is established based on the reordered ma- 
trix of FIG. 10; 

FIG. 1 2 shows a reordered matrix that is created by 
reordering elements of the coefficient matrix of FIG. 
9 in accordance with the matrix reordering process 
of FIG. 4; 

FIG. 13 shows a task graph designating a critical 
path that is established based on the reordered ma- 
trix of FIG. 12; 

FIG. 14 is a flowchart showing a conventional matrix 
reordering process in accordance with the minimum 
degree method; 

FIG. 15 shows a reordered matrix that is created by 
reordering elements of the matrix of FIG. 9 in ac- 
cordance with the conventional matrix reordering 
process of FIG. 14; 

FIG. 16 shows a task graph designating a critical 
path that is established based on the reordered ma- 
trix of FIG. 15; and 

FIG. 17 is a drawing that is used to explain upper 
and lower triangular matrices that are made by split- 
ting the matrix using a diagonal line. 
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DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

[0012] This invention will be described in further detail 
by way of examples with reference to the accompanying 
drawings. 

[0013] FIG. 1 shows an overall flow of steps for an 
electronic circuit simulation method using matrix reor- 
dering in accordance with an embodiment of the inven- 
tion. The above method is applied to an electronic circuit 
simulation system that is actualized by a computer and 
its software programs. Namely, the computer uses mul- 
tiple processors that operate in parallel to produce so- 
lutions of linear simultaneous equations in accordance 
with the software programs. The software programs are 
provided or distributed by computer-readable storage 
media or by computer networks. 
[0014] Using an input device (not shown), the system 
inputs circuit description of an electronic circuit that is a 
subject of simulation in step S201 . The electronic circuit 
contains electronic elements, each of which is ex- 
pressed by a linear equation. Hence, there are provided 
linear simultaneous equations with respect to all elec- 
tronic elements contained in the electronic circuit. In 
step S202, positions of non-zero elements selected 
from among coefficients of the linear simultaneous 
equations are registered with a matrix A. If the electronic 
circuit is expressed by linear simultaneous equations of 
N elements, the matrix A is made as 'N rows X — N col- 
umns'. In step S202, coefficients of the simultaneous 
equations are set to elements of the matrix A. In step 
S203, the system performs a matrix reordering process 
in which using matrix structure information representing 
values of elements of the matrix A, the system performs 
replacement between rows and columns in arrange- 
ment of elements of the matrix A. Herein, an order of the 
matrix A is determined to minimize a processing time of 
step S206 regarding a solving process of linear simul- 
taneous equations. In general, a combination of ele- 
ments of row i and elements of column i is called 'pivot 
i\ So, the step S203 proceeds to rearrangement of piv- 
ots contained in the matrix A. Namely, arrangement of 
elements of the matrix A is completely changed by si- 
multaneous replacement between rows and columns in 
the matrix A. 

[0015] Next, the flow proceeds to step S204 in which 
the system produces an optimum time split width A t for 
the electronic circuit simulation. So, the system sequen- 
tially determines time points by using the time split width 
A t. For example, a next time point tn+1 is determined 
using A t from a present time tn, where tn+1=tn+ A t. In 
step S203, the system provides matrix structure infor- 
mation representing new arrangement of rows and col- 
umns of the matrix A on which rearrangement of ele- 
ments is completed. Based on the matrix structure in- 
formation, non-linear characteristics of the electronic el- 
ements at a given time are linearized to create Newton 
equations (namely, linear simultaneous equations) in 



step S205. 

[0016] In step S206, the system uses Gaussian elim- 
ination processes to solve the Newton equations that 
are created in the step S205. Concretely speaking, the 

5 system operates multiple processors in parallel to solve 
the Newton equations using the Gaussian elimination 
processes. In step S207, a decision is made as to 
whether solutions to the Newton equations are con- 
verged or not. If not converged, the flow returns to step 

10 S205 in which the system re-creates Newto n equations, 
then, the flow proceeds to the solving process of step 
S206. If the solutions are converged, the flow proceeds 
to step S208 in which a decision is made as to whether 
a present time tn reaches a last time or not. If so, in other 

15 words, if tn & (last time), the system outputs calculation 
results. If the present time tn does not reach the last 
time, the system updates tn by tn+1 , then, the flow pro- 
ceeds to step S204. So, the system repeatedly performs 
a series of steps S204 to S207. 

20 [0017] The aforementioned electronic circuit simula- 
tion method is disclosed by various papers and docu- 
ments such as Japanese Patent Unexamined Publica- 
tion No. Hei 10-11475 and Japanese Patent Unexam- 
ined Publication No. Hei 7-129637, each of which is en- 

25 titled "Circuit Simulation Method", for example. 

[0018] The Gaussian elimination process of step 
S206 can be regarded as one kind of successive elimi- 
nation in which variables of simultaneous equations are 
successively eliminated one by one. Herein, a process 

30 for eliminating one variable is defined as a task. There 
exists some restriction among tasks in execution order. 
Such restriction can be represented by a digraph, which 
is called a "task graph" shown in FIG. 1 1 . For example, 
there is a restriction in which task i should be executed 

35 after completion of task j. In the task graph of FIG. 11 , 
such a restriction is translated into existence of a side 
that directs from a vertex j to a vertex i. FIG. 11 shows 
a number of vertices, each of which is represented by a 
certain integral number encompassed by a circle. In ad- 

40 dition, sides are shown by arrows, each of which repre- 
sents a direction from one vertex to another vertex. Fur- 
ther, the processing time that is required for completion 
of the task is given as a weight for the task. Such a 
weight is shown by an integral number that is described 

45 in proximity to the vertex corresponding to the task. 
[0019] Among paths each of which lies from one ver- 
tex to another in the task graph, a longest path is called 
a critical path. The parallel processing time of the Gaus- 
sian elimination process of step S206 depends upon a 

50 length of the critical path, which is a sum of weights of 
all vertices contained in the path. Therefore, the present 
embodiment performs the matrix reordering process 
(see step S203) such as to reduce the length of the crit- 
ical path as short as possible. 

55 [0020] Next, brief descriptions will be given with re- 
spect to "matrix graphs", "degrees of pivots" and "critical 
paths of pivots", which are required for explaining the 
matrix reordering process shown in FIG. 1 - The matrix 
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graph is a non-directed graph in which a side exists be- 
tween vertices i and j only when an element a (i j) that 
lies in row i and column j of the matrix A is not zero. For 
example, a matrix of FIG. 6 is translated into a matrix 
graph of FIG. 7. Generally speaking, a term "degree" 5 
denotes a number of sides connected to a vertex in a 
graph. In the case of the matrix graph, each vertex cor- 
responds to a pivot, hence, a degree of the vertex is 
called a degree of the pivot. In the matrix of FIG. 6, non- 
zero elements are denoted by marks x and O. Con- 10 
cretely speaking, numbers encompassed by circles, 
which are arranged on a diagonal line in the matrix, de- 
note numbers of pivots which are originally designated 
prior to the matrix reordering. In addition, blank positions 
represent zero elements, while a symbol 'F denotes a is 
non-zero element that is added by the Gaussian elimi- 
nation process (see step S206) with respect to each piv- 
ot. 

[0021] Suppose that during the matrix reordering of a 
square matrix of N x N in progress, pivots 1 to i are 20 
determined in order (or rearranged in position) while oth- 
er pivots i+1 to N are not determined in order. In that 
case, it is possible to perform the Gaussian elimination 
process for eliminating variables 1 to i. That is, it is pos- 
sible to create a task graph containing i tasks ranging 25 
from task 1 to task i. In the task graph, a longest path is 
selected from among paths reaching task k and is called 
a critical path of task k, which can be also called a critical 
path of pivot k. In addition, a sum of weights of all verti- 
ces contained in the critical path of task k is called a 30 
critical path length of task k or pivot k. 
[0022] FIG. 2 shows a flow of the matrix reordering 
process, which is performed in step S203 of the elec- 
tronic circuit simulation method shown in FIG. 1 . In step 
S300, a variable i is set to '1 ' while a variable N is used 35 
to designage a size of the matrix A. In step S302, initial- 
ization is performed on array variable len(k) as len(k): 
=0 with respect to all variables k (where 1 ^k^N). Here- 
in, the array variable len(k) stores a critical path length 
of pivot k that varies at each time point. In step S304, a 40 
minimum degree mindeg is detected from among de- 
grees of pivots i to N. I n step S306, pivot p whose degree 
is under mindeg+ a and whose critical path length is 
minimum is selected from among the pivots i to N. Here- 
in, a is a constant value that is zero or more. Recom- 4$ 
mended values being selected for a range from '1 ' to *3\ 
Therefore, content of the step S306 is summarized such 
that one pivot or plural pivots each having an approxi- 
mate minimum degree within a prescribed range count- 
ed from the minimum degree mindeg are examined to 50 
extract a single pivot p having a minimum (or shortest) 
critical path length. In step S308, replacement is per- 
formed between pivots p and i. 
[0023] Next, the system picks up all elements (each 
denoted by 'e') whose values are renewed by the elim- 55 
ination process of pivot i after the aforementioned re- 
placement. In step S310, a critical path length len(j) of 
pivot j to which the element e belongs are renewed as 



len(j):=max{len(i)+w(i), len(j)}. That is, critical path 
lengths of pivots that are influenced by the replacement 
(see step S308) between pivots p and i are updated by 
newest values in step S31 0. Herein, max(x,y) is a func- 
tion for selecting bigger one of variables x and y. In step 
S312, non-zero elements that are produced by the elim- 
ination process of pivot i are added to the matrix A. In 
step S31 4, the variable i is updated by adding *1 ' thereto. 
In step S31 6, a decision is made as to whether the var- 
iable i exceeds a variable N or not. Until the variable i 
exceeds the variable N, a series of steps S304 to S312 
are repeatedly performed. Thus, the matrix reordering 
process is completed with respect to all pivots contained 
in the square matrix of NxN. 

[0024] In the above, w(i) is described as a function to 
provide weight (or processing time) of pivot i. The func- 
tion w(i) can be defined as an operation to obtain a total 
number of calculations such as multiplication and divi- 
sion, which are effected in the Gaussian elimination 
process of pivot i in step S206 shown in FIG. 1 . 
[0025] Suppose that the Gaussian elimination proc- 
ess is performed in accordance with a program shown 
in FIG. 5, wherein it is possible to generalize calculation 
of w(i) as w(0=(deg(i)+1) . deg(i). That is, in an elimina- 
tion process of variable k t weight w(k) is defined as a 
result of calculations in which a degree of pivot k (i.e.; 
deg(k)) is added with T and is then multiplied by deg 
(k), namely w(k)=(deg(k)+1) • deg(k). There are various 
kinds of description methods for describing programs of 
the Gaussian elimination process. If a different program 
is described for the Gaussian elimination process, it is 
necessary to adequately change the function w(i) to re- 
flect the processing time defined by the different de- 
scription of the program. 

[0026] Next, detailed contents of the steps S304 and 
S306 shown in FIG. 2 will be described with reference 
to FIG. 3. That is, FIG. 3 shows a flow of steps corre- 
sponding to details of the steps S304 and S306, wherein 
steps S402 to S412 correspond to the step S304 while 
steps S41 4to S426 correspondto the step S306. In step 
S402, initialization is performed on a variable mindeg 
storing a minimum degree. In step S404, a variable] is 
set to a variable i. Then, a series of steps S406 to S412 
are repeatedly performed until the variable j reaches a 
variable N. Specifically, in step S406, a decision is made 
as to whether a degree of pivot j (deg(j)) is smaller than 
the variable mindeg or not. If so, the flow proceeds to 
step S408 in which the variable mindeg is renewed by 
the degree deg(j). In step S410, the variable j is in- 
creased by T, namely j:=j+1 . In step S412, a decision 
is made as to whether the increased variable j is greater 
than the variable Nor not. If not, the flow returns to step 
S406, so that a series of steps S406 to S412 are to be 
performed again. Due to the aforementioned steps of 
the step S304, a minimum degree is stored in the vari- 
able mindeg. 

[0027] In the step S306 shown in FIG. 3, initialization 
is performed on a variable minlen that stores a minimum 
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critical path length in step S414. In step S416, the var- 
iable j is set to the variable i. Then, a series of steps 
S418 to S426 are repeatedly performed until the varia- 
ble j reaches the variable N. Specifically, in step S41 8, 
a decision is made as to whether the degree of pivot j s 
(deg(j)) is smaller than mindeg+ a or not. If so, the flow 
proceeds to step S420 in which a decision is made as 
to whether a critical path length of pivot j (lenfl)) is small- 
er than the variable minlen or not. If so, the flow pro- 
ceeds to step S422 in which the variable minlen is re- 10 
newed by the critical path length len(j) t and a variable p 
is set to the variable j. After at least one of the steps 
S418 to S422 is completed, the flow proceeds to step 
S424 in which the variable j is increased by '1 ', namely 
j:=j+1 • In step S426, a decision is made as to whether is 
the variable j is greater than the variable N or not. If not, 
the flow returns to step S418, so that a series of steps 
S418 to S426 are to be performed again. Due to the 
aforementioned steps of the step S306, the variable p 
stores a number of a pivot whose degree is under min- 20 
deg+ a and whose critical path length is minimum. In 
addition, the variable minlen stores the minimum critical 
path length. 

[0028] Next, a concrete example of the matrix reor- 
dering process, which is shown in step S203 of FIG. 1 25 
and details of which are shown in FIG. 2, will be de- 
scribed with reference to Figures 8 to 11 . FIG. 8 shows 
an example of an electronic circuit to be simulated. FIG. 
9 shows an example of a coefficient matrix representing 
simultaneous equations corresponding to the electronic 30 
circuit. FIG. 1 0 shows an example of a reordered matrix 
that is created by effecting matrix reordering on the co- 
efficient matrix of FIG. 9. FIG. 11 shows a task graph 
representing the reordered matrix of FIG. 1 0. 
[0029] That is, FIG. 8 shows an electronic circuit that 35 
is configured by nine electronic elements D to D9, which 
are connected together via seven nodes N1 toN7. Here- 
in, voltage variations of the nodes N1 to N7 are repre- 
sented by variables, by which linear simultaneous equa- 
tions are created. FIG. 9 shows a coefficient matrix that 40 
represents the linear simultaneous equations with re- 
spect to the seven variables. The coefficient matrix con- 
tains a certain number of pivots, wherein pivot i is con- 
figured to cope with a node number i. In the coefficient 
matrix of FIG. 9, non-zero elements are designated by 4s 
marks x, and zero elements are located at blank posi- 
tions. For convenience' sake, numbers of the elements 
are omitted in the coefficient matrix of FIG. 9. In addition, 
marks O also designate non-zero elements, wherein in- 
tegral numbers encompassed by circles designate n urn- so 
bers of pivots before matrix reordering. Weight w(i) for 
a Gaussian elimination process of a variable i (namely, 
task i) is produced by multiplication between '(a degree 
of pivot i)+V and 'a degree of pivot i', namely w(i)=(deg 
(i)+1) • deg(i). The following description is made under ss 
a prescribed condition where a=0. 
[0030] In step S300 of FIG. 2, a variable i is set to '1 \ 
Herein, pivot i is determined as a pivot whose degree is 



minimum (where o=0). Then, the system searches a 
pivot whose critical path length is minimum. In FIG. 9, 
all of pivots 1 to 7 have a minimum degree '1 ' and a min- 
imum critical path length '2'. Among those pivots 1-7, 
the system chooses pivot 1 whose number is the small- 
est. Hence, a variable p is set to '1' in step S306. Be- 
cause of i=p, the system does not actually perform re- 
placement of pivots in step S308. Even if a Gaussian 
elimination process is performed on the pivot 1 , non-ze- 
ro elements are not newly produced (see step S312). 
[0031] If i=3, pivot 3 having a minimum degree '2' and 
a minimum critical path length '6* is selected, so that the 
variable p is set to '3'. By elimination of the pivot 3, non- 
zero elements a(4,7) and a(7,4) are newly produced and 
are added to the matrix in step S31 2. 
[0032] If i=4, pivot 5 having a minimum degree '2' and 
a minimum critical path length '6' is selected, so that the 
variable p is set to '5'. By elimination of the pivots, non- 
zero elements a(4,6) and a(6,4) are newly produced and 
are added to the matrix in step S312. 
[0033] If i=5, pivot 4 having a minimum degree '2' and 
a minimum critical path length '12' (=6+6) is selected, 
so that the variable p is set to '4\ By elimination of the 
pivot 4, non-zero elements are not newly produced (see 
step S312). 

[0034] If i=6, pivot 6 having a minimum degree '1 ' and 
a minimum critical path length '14' (=6+6+2) is selected, 
so that the variable p is set to '6\ By elimination of the 
pivot 6, non-zero elements are not newly produced (see 
stepS312). 

[0035] Due to the aforementioned processes, the co- 
efficient matrix of FIG. 9 is subjected to matrix reordering 
and is transformed to a reordered matrix shown in FIG. 
10. A task graph for the reordered matrix is shown in 
FIG. 11. Herein, a critical path is established among 
tasks 3, 4, 6 and 2, so that task 2 having a degree '0' is 
the last task that is executed last in order. In addition, 
FIG. 11 shows that a critical path length is '14', and a 
sum of weights is '24'. 

[0036] If the conventional matrix reordering shown in 
FIG. 14 is applied to the coefficient matrix of FIG. 9, it is 
possible to create a reordered matrix shown in FIG. 15. 
FIG. 16 shows a task graph for the reordered matrix of 
FIG. 15 to establish a critical path among tasks 3, 4, 5, 
6 and 2, so that task 2 is the last task that is executed 
last in order. In addition, FIG. 16 shows that a critical 
path length is '20', and a sum of weights is '24'. 
[0037] As described above, the matrix reordering 
process of step S203 of FIG. 1 is described with refer- 
ence to Figures 2 and 3, wherein the matrix reordering 
process is performed in consideration of the approxi- 
mate least degree (=minimum degree + a) and critical 
path length. As compared with the conventional tech- 
nique, the present embodiment can considerably re- 
duce a total processing time in parallel processing using 
multiple processors that operate in parallel with respect 
to tasks. 

[0038] Next, a description will be given with respect 
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to the matrix reordering process (step S203) in accord- 
ance with another embodiment with reference to FIG. 4. 
FIG. 4 shows a flow of steps of the matrix reordering 
process that is contained in the aforementioned elec- 
tronic circuit simulation method shown in FIG. 1. Rough- 5 
ly speaking, the matrix reordering process of FIG. 4 is 
made by a combination of two reordering methods, each 
of which can be selectively activated. That is, a first re- 
ordering method is an approximate minimum degree 
and minimum critical path method (see step S8) that is 10 
similar to the aforementioned method shown in FIG. 2 
in which reordering is performed in consideration of the 
approximate minimum degree and minimum critical 
path length, and a second reordering method corre- 
sponds to one of the conventional matrix reordering 15 
methods, namely a nested dissection method (see step 
S10). 

[0039] In FIG. 4, step S8 representing the approxi- 
mate minimum degree and minimum critical path meth- 
od is configured by seven steps S1 to S7. In step S1 , 20 
initialization is performed so that a variable i is set to "1 '. 
In step S2, the system selects from among pivots that 
are not determined in order (or arrangement) in a given 
coefficient matrix A, a pivot (whose number is p) that 
has an approximate minimum degree (which is under 25 
minimum degree + a, where a recommended value for 
a ranges from '1 ' to '3') and a minimum critical path. In 
step S3, replacement is performed between pivots p and 
i, so that new pivot i is to be determined. In step S4, non- 
zero elements, which are produced by a Gaussian elim- 30 
ination process of pivot i, are added to the coefficient 
matrix A. A series of the aforementioned steps are re- 
peatedly performed until a degree of pivot i (deg(i)) ex- 
ceeds a constant p in step S5, or they are repeatedly 
performed until all pivots are determined in order (or ar- 35 
rangement) in step S7. 

[0040] As compared with the nested dissection meth- 
od, the approximate minimum degree and minimum crit- 
ical path method of step S8 is advantageous because it 
is possible to further reduce a critical path length when 40 
an average degree of pivots contained in the given ma- 
trix is approximately under '3'. However, if the average 
degree becomes greater than '3', the approximate min- 
imum degree and minimum critical path method tends 
to produce a longer critical path length as compared with 45 
the nested dissection method. For this reason, during 
execution of the approximate minimum degree and min- 
imum critical path method in progress, when a degree 
of the selected pivot exceeds the constant (J (where 0 is 
a natural number), the system automatically changes so 
from the approximate minimum degree and minimum 
critical path method to the nested dissection method. In- 
cidentally, it is preferable to set a range between '3' and 
'10' for the constant p. 

[0041] If all pivots are not yet determined in order al- 55 
though the matrix reordering is effected on the coeffi- 
cient matrix A in accordance with the approximate min- 
imum degree and minimum critical path method of step 



S8, the flow proceeds to step S9 in which a matrix graph 
is created for a partial matrix R with respect to non-elim- 
inated pivots of the matrix A. In step S10, the matrix 
graph is used to reorder elements of the partial matrix 
R in accordance with the nested dissection method. In 
the nested dissection method, the matrix graph is recur- 
sively and repeatedly split into two sections, wherein re- 
ordering is performed such that a pivot corresponding 
to a vertex on a split boundary comes last. Until all pivots 
are reordered, splitting the matrix graph is recursively 
and repeatedly performed. In step S11, non-zero ele- 
ments, which are produced by the Gaussian elimination 
process of the partial matrix R, are added to the coeffi- 
cient matrix A in step S11. Thus, it is possible to com- 
plete the reordering process of the coefficient matrix A. 
[0042] Next, a concrete example of the matrix reor- 
dering process, which is shown in step S203 of FIG. 1 
and details of which are shown in FIG. 4, will be de- 
scribed with reference to Figures 6 to 9 and Figures 12, 
13. As similar to the foregoing embodiment, the matrix 
reordering process is performed using the electronic cir- 
cuit of FIG. 8 that is configured by nine electronic ele- 
ments D1 -D9 via seven nodes N1 -N7. That is, the matrix 
reordering process is performed on a coefficient matrix 
shown in FIG. 9 that is created based on linear simulta- 
neous equations, which express seven variables corre- 
sponding to the seven nodes N1-N7 of the electronic 
circuit. In addition, weight w(i) for task i (i.e., Gaussian 
elimination process of variable i) is produced by multi- 
plication between '(a degree of pivot i)+1 ' and 'a degree 
of pivot i', namely w(i)=(deg(i)+1) - deg(i). The following 
description is made under the condition where a =0, p 
=1. 

[0043] In the approximate minimum degree and min- 
imum critical path method of step S8 shown in FIG. 4, 
a flow firstly proceeds to step S1 in which a variable i is 
set to '1'. Herein, pivot i is determined as a pivot whose 
degree is minimum (where a =0). Then, the system 
searches a pivot whose critical path length is minimum 
with respecttothe pivot i. In the coefficient matrix of FIG. 
9, both of pivots 1 , 7 have a minimum degree '1 ' and a 
minimum critical path length '2'. Within those pivots 1 , 
7, the system selects pivot 1 whose number is smaller, 
so that a variable p is set to '1 ' in step S2. Because of 
i=p, replacement is not actually performed between piv- 
ots p and 

i. In addition, non-zero elements are not produced by a 
Gaussian elimination process of pivot 1 (see step S4). 
[0044] If fc=2, the aforementioned processes are re- 
peated on elements of the matrix that are not deter- 
mined in order. Herein, pivot 7 has a minimum degree 
and a minimum critical path length, so the variable p is 
set to 7*. In step S3, replacement is performed between 
pivots 2 and 7. However, non-zero elements are not pro- 
duced by a Gaussian elimination process of pivot 7 (see 
step S4). 

[0045] If i=3, pivot 3 has a minimum degree '2' and a 
minimum critical path length '6', so the variable p is set 
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to '3'. By elimination of pivot 3, non-zero elements a(4,7) 
and a(7,4) are newly produced and are added to the ma- 
trix in step S4. Since the degree of pivot 3 is '2' that is 
greater than p =1 , the matrix reordering process based 
on the approximate minimum degree and minimum crit- 
ical path method is ended, so the flow proceeds from 
step S5 to step S9. 

[0046] In the above, a partial matrix R as shown in 
FIG. 6 is made by pivots of the matrix A that are not 
determined in order. A matrix graph shown in FIG. 7 is 
created for the partial matrix R of FIG. 6 in step S9. The 
nested dissection method splits the matrix graph into 
two sections, so that the pivots of the partial matrix R 
are determined in order such that a pivot corresponding 
to a node on a split boundary comes last. Splitting the 
matrix graph is recursively and repeatedly performed 
until all pivots are completely determined in order. In the 
case of the partial matrix R, it is sufficient to split the 
matrix graph one time. That is, pivots 4 and 6 corre- 
spond to vertices on a split boundary P1 , hence, pivots 
remained in the partial matrix R are reordered as pivots 
5, 2, 4 and 6. FIG. 12 shows a reordered matrix that is 
a final result of the matrix reordering process described 
above. FIG. 13 shows a task graph for use in the Gaus- 
sian elimination process of the reordered matrix of FIG. 
12. Numbers of vertices included in the task graph cor- 
respond to numbers of pivots included in the reordered 
matrix, wherein a critical path is established among piv- 
ots 3, 2, 4 and 6. In the task graph, a critical path length 
is '14'. This proves that the matrix reordering process of 
FIG. 4 can reduce the critical path length as compared 
with the conventional matrix reordering process of FIG. 
14. 

[0047] The aforementioned matrix reordering proc- 
esses of Figures 3 and 4 (i.e., step S203 shown in FIG. 
1) are described in connection with symmetric matrices. 
Of course, it is possible to apply the matrix reordering 
processes to non-symmetric matrices. In that case, it is 
necessary to introduce three points of modification to 
the matrix reordering processes of the foregoing em- 
bodiments, which will be described below. 
[0048] A first point of modification is commonly ap- 
plied to both of the matrix reordering processes of Fig- 
ures 2 and 4, namely a necessary condition for selecting 
a pivot having a minimum critical path length. The afore- 
mentioned matrix reordering processes use a condition 
where a degree deg(i) is under the prescribed value a, 
which is modified to another condition where a param- 
eter *jNZLC(i) NZUR(i) is under the prescribed value 
a. FIG. 17 shows that a matrix is split using a diagonal 
line into two triangular matrices L and U, both of which 
do not contain elements on the diagonal line. In the 
aforementioned parameter, NZLC(i) represents a 
number of non-zero elements of column i within the low- 
er triangular matrix L, and NZUR(i) represents a number 
of non-zero elements of row i within the upper triangular 
matrix U. 

[0049] Second and third points of modification are ex- 



clusively applied to the matrix reordering process of FIG. 
4 containing a combination of the approximate minimum 
degree and minimum critical path method and the nest- 
ed dissection method, which are selectively used in re- 
5 sponse to the prescribed switch condition (see step S5) . 
The approximate minimum degree and minimum critical 
path method shown in FIG. 2 can be applied to non-sym- 
metric matrices by merely introducing the first point of 
modification. The second point of modification is a 

10 change of the switch condition forswitching between the 
approximate minimum degree and minimum critical 
path method of step S8 and the nested dissection meth- 
od of step S10. That is, the switch condition of step S5 
where deg(i)>p is modified using the aforementioned 

* 5 p arameter to another condition where 
JVZLC(i) NZUR(i) > p. The third point of modification is 
a change of application of the nested dissection method 
of step St 0. That is, the nested dissection method is not 
applied directly to a matrix graph of the non-symmetric 

20 partial matrix R but is applied to a matrix graph of the 
symmetric matrix (R+RT) that is created by transposi- 
tion of the non-symmetric matrix R. Herein, RT repre- 
sents a transpose matrix that is created by interchang- 
ing rows and columns of the matrix R. By introducing 

25 the aforementioned second and third points of modifi- 
cation, it is possible to apply the matrix reordering proc- 
ess of FIG. 4 to non-symmetric matrices. 
[0050] In the aforementioned modification that is ef- 
fected on the matrix reordering process of FIG. 4, prior 

30 to reordering of elements of a non-symmetric matrix A, 
it is necessary to create a symmetric matrix (A+AT) by 
the non-symmetric matrix A and its transpose matrix AT. 
Then, elements of the symmetric matrix (A+AT) are re- 
ordered by the approximate minimum degree and min- 

35 imum critical path method of step S8 and the nested dis- 
section method of step S1 0. 

[0051 ] Application of the matrix reordering processes 
of this invention is not necessarily limited to the elec- 
tronic circuit simulation method. In other words, this in- 

40 vention can be applied to any other types of data 
processing methods including solving processes for 
solving linear simultaneous equations by the parallel 
processing of Gaussian elimination processes. 
[0052] In short, this invention is improved as com- 

45 pared with the conventional methods because of con- 
siderable reduction of processing times for the parallel 
processing of Gaussian processes in which critical path 
lengths are sufficiently reduced in task graphs. 
[0053] As this invention may be embodied in several 

50 forms without departing from the spirit of essential char- 
acteristics thereof, the present embodiments are there- 
fore illustrative and not restrictive, since the scope of the 
invention is defined by the appended claims rather than 
by the description preceding them, and all changes that 

55 fall within metes and bounds of the claims, or equiva- 
lence of such metes and bounds are therefore intended 
to be embraced by the claims. 
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Claims 

1 . A matrix reordering method for reordering elements 
of a coefficient matrix whose structure corresponds 

to coefficients of linear simultaneous equations $ 
whose solutions are to be produced by parallel 
processing of processors of a computer in accord- 
ance with Gaussian elimination, said matrix reor- 
dering method comprising the steps of: 

10 

based on a number of non-zero elements in- 
cluded in the coefficient matrix and an accumu- 
lative processing time of the Gaussian elimina- 
tion of the coefficient matrix, determining a first 
combination of row and column and a second 15 
combination of row and column, which are se- 
lected from among rows and columns of the co- 
efficient matrix; and 

performing replacement of elements between 
the first combination of row and column and the 20 
second combination of row and column within 
the coefficient matrix. 

2. A matrix reorderi ng method for reordering elements 

of a coefficient matrix whose structure corresponds 25 
to coefficients of linear simultaneous equations 
whose solutions are to be produced by parallel 
processing of processors of a computer in accord- 
ance with Gaussian elimination, said matrix reor- 
dering method comprising the steps of: 30 

based on a number of non-zero elements in- 
cluded in the coefficient matrix and lengths of 
critical paths created by the Gaussian elimina- 
tion of the coefficient matrix, determining a first 35 
combination of row and column and a second 
combination of row and column, which are se- 
lected from among rows and columns of the co- 
efficient matrix; and 

performing replacement of elements between *o 
the first combination of row and column and the 
second combination of row and column within 
the coefficient matrix. 

3. A matrix reordering method according to claim 1 or 
2 further comprising the step of: 

in accordance with a prescribed condition, se- 
lectively performing either the replacement of ele- 
ments between the first and second combinations 
of rows and columns or secondary replacement of so 
elements between a third combination of row and 
column and a fourth combination of row and col- 
umn , which are selected based on symmetry of the 
coefficient matrix. 

55 

4. A matrix reordering method according to claim 3 fur- 
ther comprising the step of: 

creating a symmetric coefficient matrix by 



transposition of a non-symmetric coefficient matrix 
that is given as the coefficient matrix, so that the 
secondary replacement is performed on elements 
of the symmetric coefficient matrix between the third 
and fourth combinations of rows and columns, 
which are selected based on symmetry of the sym- 
metric coefficient matrix. 

5. A matrix reordering apparatus for reordering ele- 
ments of a coefficient matrix whose structure corre- 
sponds to coefficients of linear simultaneous equa- 
tions whose solutions are to be produced by parallel 
processing of processors of a computer in accord- 
ance with Gaussian elimination, said matrix reor- 
dering apparatus comprising: 

a determinator for based on a number of non- 
zero elements included in the coefficient matrix 
and an accumulative processing time of the 
Gaussian elimination of the coefficient matrix, 
determining a first combination of row and col- 
umn and a second combination of row and col- 
umn, which are selected from among rows and 
columns of the coefficient matrix; and 
a replacer for performing replacement of ele- 
ments between the first combination of row and 
column and the second combination of row and 
column within the coefficient matrix. 

6. A matrix reordering apparatus for reordering ele- 
ments of a coefficient matrix whose structure corre- 
sponds to coefficients of linear simultaneous equa- 
tions whose solutions are to be produced by parallel 
processing of processors of a computer in accord- 
ance with Gaussian elimination, said matrix reor- 
dering apparatus comprising: 

a determinator for based on a number of non- 
zero elements included in the coefficient matrix 
and lengths of critical paths created by the 
Gaussian elimination of the coefficient matrix, 
determining a first combination of row and col- 
umn and a second combination of row and col- 
umn, which are selected from among rows and 
columns of the coefficient matrix; and 
a replacer for performing replacement of ele- 
ments between the first combination of row and 
column and the second combination of row and 
column within the coefficient matrix. 

7. An electronic circuit simulation method using matrix 
reordering for reordering elements of a coefficient 
matrix that represents electronic elements of a giv- 
en electronic circuit by linear simultaneous equa- 
tions whose solutions are to be produced by parallel 
processing of processors of a computer in accord- 
ance with Gaussian elimination, said electronic cir- 
cuit simulation method comprising the steps of: 
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based on a number of non-zero elements in- 
cluded in the coefficient matrix and an accumu- 
lative processing time of the Gaussian elimina- 
tion of the coefficient matrix, determining a first 
combination of row and column and a second s 
combination of row and column, which are se- 
lected from among rows and columns of the co- 
efficient matrix; and 

performing replacement of elements between 
the first combination of row and column and the io 
second combination of row and column within 
the coefficient matrix. 

8. An electronic circuit simulation method using matrix 
reordering for reordering elements of a coefficient '5 
matrix that represents electronic elements of a giv- 
en electronic circuit by linear simultaneous equa- 
tions whose solutions are to be produced by parallel 
processing of processors of a computer in accord- 
ance with Gaussian elimination, said electronic cir- 20 
cuit simulation method comprising the steps of: 

based on a number of non-zero elements in- 
cluded in the coefficient matrix and lengths of 
critical paths created by the Gaussian elimina- 25 
tion of the coefficient matrix, determining a first 
combination of row and column and a second 
combination of row and column, which are se- 
lected from among rows and columns of the co- 
efficient matrix; and 30 
performing replacement of elements between 
the first combination of row and column and the 
second combination of row and column within 
the coefficient matrix. 

35 

9. An electronic circuit simulation method according to 
claim 7 or 8 further comprising the step of: 

in accordance with a prescribed condition, se- 
lectively performing either the replacement of ele- 
ments between the first and second combinations 40 
of rows and columns or secondary replacement of 
elements between a third combination of row and 
column and a fourth combination of row and col- 
umn, which are selected based on symmetry of the 
coefficient matrix. 45 

1 0. An electronic circuit simulation method according to 
claim 9 further comprising the step of: 

creating a symmetric coefficient matrix by 
transposition of a non-symmetric coefficient matrix so 
that is given as the coefficient matrix, so that the 
secondary replacement is performed on elements 
of the symmetric coefficient matrix between the third 
and fourth combinations of rows and columns, 
which are selected based on symmetry of the sym- 55 
metric coefficient matrix. 

1 1 . An electronic circuit simulation apparatus using ma- 



trix reordering for reordering elements of a coeffi- 
cient matrix that represents electronic elements of 
a given electronic circuit by linear simultaneous 
equations whose solutions are to be produced by 
parallel processing of processors of a computer in 
accordance with Gaussian elimination, said elec- 
tronic circuit simulation apparatus comprising: 

a determinator for based on a number of non- 
zero elements included in the coefficient matrix 
and an accumulative processing time of the 
Gaussian elimination of the coefficient matrix, 
determining a first combination of row and col- 
umn and a second combination of row and col- 
umn, which are selected from among rows and 
columns of the coefficient matrix; and 
a replacer for performing replacement of ele- 
ments between the first combination of row and 
column and the second combination of row and 
column within the coefficient matrix. 

12. An electronic circuit simulation apparatus using ma- 
trix reordering for reordering elements of a coeffi- 
cient matrix that represents electronic elements of 
a given electronic circuit by linear simultaneous 
equations whose solutions are to be produced by 
parallel processing of processors of a computer in 
accordance with Gaussian elimination, said elec- 
tronic circuit simulation apparatus comprising the 
steps of: 

a determinator for based on a number of non- 
zero elements included in the coefficient matrix 
and lengths of critical paths created by the 
Gaussian elimination of the coefficient matrix, 
determining a first combination of row and col- 
umn and a second combination of row and col- 
umn, which are selected from among rows and 
columns of the coefficient matrix; and 
a replacer for performing replacement of ele- 
ments between the first combination of row and 
column and the second combination of row and 
column within the coefficient matrix. 

13. A matrix reordering method for reordering elements 
of a coefficient matrix created based on coefficients 
of linear simultaneous equations whose solutions 
are to be produced by parallel processing of proc- 
essors of a computer in accordance with Gaussian 
elimination, said matrix reordering method compris- 
ing the steps of: 

selecting from among pivots included in the co- 
efficient matrix a first pivot whose degree cor- 
responding to a number of non-zero elements 
is under a threshold; 

selecting from among the pivots included in the 
coefficient matrix a second pivot whose critical 
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path length is minimum; 
performing replacement of elements between 
the first pivot and the second pivot within the 
coefficient matrix; and 

adding new non-zero elements, which are new- 5 
ly produced by the Gaussian elimination of the 
first pivot, to the coefficient matrix. 

14. A matrix reordering method according to claim 13 
further comprising the step of: io 

performing reordering of a partial matrix 
whose elements are not eliminated and are select- 
ed from among the elements of the coefficient ma- 
trix in accordance with a nested dissection method, 
so that non-zero elements, which are newly pro- is 
duced by the Gaussian elimination of the partial ma- 
trix, are added to the coefficient matrix. 

15. A matrix reordering method according to claim 14 
further comprising the step of: 20 

creating a symmetric coefficient matrix by 
transposition of a non-symmetric coefficient matrix 
that is given as the coefficient matrix, so that the 
reordering is performed on the symmetric coeffi- 
cient matrix. 25 

16. A matrix reordering method according to claim 14 
wherein the reordering is started if a degree or a 
parameter of the first pivot is under a threshold. 

30 
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FIG. 1 
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FIG. 4 
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